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Abstract 

B plv decay is analyzed in the effective tfieory of heavy quark with 
infinite mass hmit. The matrix element relevant to the heavy to light vector 
meson semileptonic decays is parametrized by a set of four heavy flavor-spin 
independent universal wave functions at the leading order of effective theory. 
The form factors are calculated at the leading 1 / mq order using the light cone 
sum rule method in the framework of effective theory. \Vub\ is then extracted 
via B plv decay mode. 

PACS numbers: 13.20.-v, ll.30.Hv, 11.55.Hx, 12.15.Hh, 12.39.Hg 
Keywords: heavy to light vector transition, |T4fe|; effective theory, light cone 
sum rule 
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Typeset using REVT^ 



I. INTRODUCTION 



B — > 7r(p)Zz/ decays have been studied in the full QCD via sum rule approach 
or using appropriate models In Refs. |jlO[ and [0, the heavy to light pseudoscalar 



and vector meson decay matrix elements in the heavy quark effective theory have been 
formulated respectively up to the order of l/mq. In Ref. [0, the leading order wave 
functions of i? — *• ttIi' decay have been calculated in the effective theory of heavy quark by 
using the light cone sum rule method, and the important CKM matrix element \Vub\ has 
been extracted. 

This short paper is parallel to Ref. ||T^ and focuses on the calculation of the leading 
order wave functions of B ^ plv decay and the extraction of \Vub\ by using the light cone 
sum rules in the effective theory of heavy quark. It is organized as follows: In section 2, 
we present the analytic formulae for B — >■ plv decay in the effective theory of heavy quark. 
In section 3, the relevant light cone sum rules are derived within the framework of effective 
theory. In section 4 we present the numerical results including the extracted value of \Vub\- 
The value of \Vub\ is compared with that extracted via B ttIu decay. A short summary is 
drawed in the last section. 



II. MATRIX ELEMENT 

The B plu decay matrix element is generally written in terms of four form factors as 

< p(p, e*)|M7^(l - 7')&|5(P + q) >= -^{ruB + mp)Ai(g')e*^ + ^^M^(e* ■ (p + g))(2p + qY 

rriB + rrip 

■ (p + q))q^ - ^^^^ e^^^^^e^p + q)pp,, (2.1) 

where p and e* are the momentum and polarization vectors of the p meson, and q is the 
momentum carried by the lepton pair. 

In the effective theory of heavy quark, heavy quark expansion (HQE) can be applied to 
the matrix element p!0|- pr3| , p!6| , p!7| . The normalization of the matrix elements in full QCD 



and in the effective theory is [l^JT3|JT^JT7|l 



< p{p,e*)\uVh\B >= -^{< p{p,e*)\uTQ+\B, > +0(1^)}, (2.2) 

^Ab 

where Ab = rriB — mb, and 

A = lim Ab (2.3) 

mQ—>oo 

is the heavy flavor independent binding energy reflecting the effects of the light degrees of 
freedom in the heavy hadron. is the effective heavy quark field. 

Based on the heavy quark symmetry (HQS), one can parametrize the leading order 



matrix element in the effective theory of heavy quark on the rhs. of eq. (|2.2|) as |11 



< p{p,e*)\uTQt\B, >= -iTT[n{v,p)TM,] (2.4) 
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with the decomposition of Q{v,p): 



niv,p) = L,iv p)^* + L2iv- p){v ■ e*) + [L,iv ■ p)^* + L,iv ■ p){v e*)]^, 



V ■ p 



(2.5) 



Aiy is the well known spin wave function associated with the heavy meson states, 



_ + ?^ f ~75 for pseudoscalar meson 
^ 2 1 / for vector meson with j 



polarization vector e^. 



(2.6) 



Li{i = 1, 2, 3, 4) are the leading order wave functions in the effective theory. Generally, they 
do not depend on the heavy quark mass but are functions of the variable v-p and the energy 
scale n as well. For the sake of simplicity we do not write down their /i dependence explicitly 
until the numerical analysis in section 
Eqs. (|0|)-(|2;|) lead to 



rriBA 
mB + rrip'^ Ab 

2{mB + m 



{L^{v p) + Lsiv p)} + ■ ■ ■ ; 



\ Ab 



2mi 



2mB{v ■ p) 



Az{q^) = 2{mB + m 



TUB A L2{v-p) _ Lsjv-p) - L^jv-p) ^ 



\ Ab 



2mi 



2mB{v ■ p) 



itlbA niB + rrip 



\ Ab mB{v-p) 



Lz{v-p) + 



(2.7) 



The dots denote higher order l/mq contributions which will not be taken into account in 
the following calculations. 



III. LIGHT CONE SUM RULE 

The form factors Ai,A2,A3 and V ( and/or Li{i = 1,2,3,4) ) are difficult to be pre- 
dicted straightforwardly from QCD calculation due to their nonperturbative nature. QCD 
sum rules, quark models and lattice simulations are the main approaches to evaluate the 
nonperturbative contributions. For heavy to light decays, the light cone sum rule approach 
is found to be more reliable than the sum rules based on short distance operator product 
expansion (OPE) and vacuum condensate input. In the light cone sum rule calculation, 
the relevant correlation functions are expanded near the light cone, and the nonperturbative 
contributions are introduced into the treatment through the light cone distribution functions 
of the mesons. 

We now consider the vacuum-p correlation function 

F^'{p,q) =i J d^xe-'P^-"" < p{p,e*)\T{u{0)Yil - l^)b{0),b{x)^l^d{x)}\0 > . (3.1) 
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Here the B meson has momentum Pb = P + Q, whereas p and q are momenta carried by 
the rho-meson and leptons. Phenomenologically, we insert a complete set of states with B 
meson quantum numbers in 

< p{p,e*)\u-f^'{l--f^)b\B >< B\bi-f^d\0 > 



F''{p,q)phe 
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m%-{p + q) 
< p{p, e*)|n7ni - l'')^H >< H\U^^d\0 > 



Due to eqs. (|2.2|) , ( |2.4|) and the B decay constant definition 



< 0\qTQ+\B, >= ^Tt[TM,], (3.3) 
at the leading order of HQE, q)phen contributes 

""^Z - -{(Li + - L2V^{e* ■ v) - (L3 - L^)p^^-—^ 

mhAB 2Ab — 2v ■ k v ■ p 

+i^e^"''"^elpavs}+ r ds ^^^ " ^' ^} + subtractions, (3.4) 
V ■ p J so s — 2v ■ k 

with /c'^ being the heavy hadron's residual momentum, k^ = — nibV^. The integral 
represents the higher resonance contributions. 

The correlator can also be calculated in theory and written as 

ds^^ — P) s) theory _^ subtractions. (3.5) 
s — 2v ■ k 

Furthermore the quark-hadron duality assumes the equality of p{v-p, s)theory and the physical 
spectral density p{y ■ p,s). As a result, equating (p.4[) and (|3.5|) yields 



00 



rrihAB 2Ab — 2v ■ k v-p v-p 

ds ^^^ p, g) subtractions. (3.6) 
s — 2v ■ k 

We now calculate the correlator ( p.l|) in the effective theory of heavy quark. When 
neglecting higher order l/mg corrections eq. ( p.l| ) can be written as 

F^ip^q) = I f d4a,e-iPBX+im,..x. ^ p(p, e*)|T{n(0)7^(l - 75)g+(0), Q+(a;)z7^rf(x)}|0 > . (3.7) 



This could then be expanded into a series in powers of the twist of light cone distribution 
functions. The p distribution functions are defined by the following matrix elements p|,p^ 



< p{p,e*)\u{0)a^^d{x)\0 > = -ifpie^pu - tlp^) j due 

J 



I {u\ 



e ■ X 



<p{p,e*)\u{Oh^d{x)\0> =fpmpp^—^ f due'^^-^(t>\\{u) 

+ fpmp{e;-p/—^) [' due'^P-^g'f\u), 
p ■ X Jo 

< p{p,e*)\u{0)j,j5d{x)\0 > = -Jpmpe^.^^e^p'^x^ f^due'^''-^g^l\u). (3.8) 



0_L and 0|j are twist 2 distribution functions of transversely and longitudinally polarized p 
mesons respectively while g'f^ and g^f* are associated with both twist 2 and twist 3 opera- 
tors. Higher twist components still lack for more precise analysis, and they are beyond our 
consideration in this paper. 

Contracting the effective heavy quark fields into the propagator dtS{—x — vt), 

eq.(|3.7|) turns into 

F^iy,uj) = - dt due-e-^-y'{f,m^p^ <P\\{u) + f^mp{e*^-p'^ )gt\u) 

2 Jo Jo p ■ V p ■ V 

+ if^m^e^'^^Pelp^vpgt\u) + f^[{v ■ p)e*^ - (v ■ e*)^ + le^^"^ elp^vp\(^^{u)} (3.9) 

with y = V ■ p and u = 2v ■ k. 

Performing a wick rotation of the t axis and using the feature of Borel transformation: 

4-)e^- = 5(A-i), (3.10) 

we get from (p.9| ) 

BPF^{y,u) = due-^^y{'^p>^[f,mMu) - g'f\u)) - f;t{vp)M^^)] 

+e*'^[f^mpg'f\u) + f^{v ■ p)Mu)] + e'^^"^eX^/5[^/,m,^?P(«) + ^f^Mu)]}■ (3.11) 



Following the approach in [|1^,|1^, we can now carry out continuous double Borel trans- 
formations on the correlator itself to produce the spectral function p{y, s). 

p{y,s) = B[j^^^^B^-^F'^iy,u;) 
zy V ■ p 

+f;^{v ■ p)M^)] + e>^''^^elp.vp[^f,m,{^/l\u)) + */>x(«)]}.=^. (3.12) 



■Ay-" ''du' 



Eq. ( p.l2| ) can be easily derived from ( p.ll| ) by first writing ^ as a derivative of the exponent 
in ( |3.11|) over u, and then using the method of integration by parts over u. 
We then get from eqs. ( |3.6| ) and ( |3.12| ) 



Uy) = 1 e^^-/^ /" rf.e-/^i/,m,[0||(«) - (..) - h^g^r\u))U^, (3.13) 

and L2{y) equals zero in the present approximation since no twist 2 distribution functions 
contribute to it. For this reason, the form factors A2 and ^3 have the same absolute value 
but opposite signs at the order considered in this paper. 
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IV. NUMERICAL RESULTS 



The functions and 0|| give the leading twist distributions in the fraction of total 
monmentum carried by the quark in transversely and longitudinally polarized mesons, re- 
spectively. They have a non-trivial scale dependence which can be described by the renor- 
malization group method |^. These distribution functions can be expanded in Gegenbauer 
polynomials C'^/'^{x) whose coefficients are renormalized multiplicatively. Namely, writting 
their scale dependence explicitly, we have 



[u.^i) = Qu{l-u)\l+ a^(ll)(/i)C^/2(2M-l)], 

n=2,4,--- 



(4.1) 



where /3o = 11 — (2/3)nj, and the one loop anomalous dimensions are p0| , |21 



8 2 n+i ^ 

^11 = -fl 1 -I- 4 V -) 

O 71+1 1 

7n =tt(1 + 4E-)- (4-2) 

^ j=2 J 

The coefficients and a| themselves are nonperturbative parameters, and have been calcu- 
lated using sum rule methods in Ref . . In the following discussions we will use the values 

i 

a^(lGeV) = 0.2±0.1, 4(lGeV) = 0.18 ± 0.10 (4.3) 

and a^^^^^ = for n ^ 2. 

The functions g^^^ and g^^^ describe transverse polarizations of quarks in the longitudi- 
nally polarized mesons. They receive contribtuions of both twist 2 and twist 3. And the twist 
2 contributions are related to the longitudinal distribtuion 0|| [u, fi) by Wandzura-Wilczek 
type relations 



2 Jo 1 — V Ju V 
(«,/^) = 2[(l-«) rdv^^ + uf'dv^^]. (4.4) 



For the energy scale fi to be used in the sum rules ( |3.13| ), we use 



fib ~ \Jm% -ml^ 2.4GeV, (4.5) 

which is an appropriate choice of scale set by the typical virtuality of the beautiful quark 
i- 

The values of the hadron quantities fp, f^, A^, A and F are needed in order to perform 
the sum rule numerical analysis. The decay constant fp has been measured in experiments 



22,23|. is the tensor coupling defined by 
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< 0\u(T^^d\p'^{p, e) >= i{epPu - euPu)fp (4.6) 

and its value has been calculated in Ref. ||^. A^, A and F are associated with heavy mesons 
and are parameters in the effective theory of heavy quark. Their values have been estimated 
consistently in Ref. |]T6[ by sum rules in the framework of effective theory. We use for these 
parameters the following values, 

fp± = (195 ± 7)MeV, f^o = (216 ± 5)MeV, = (160 ± 10)MeV, 

Ab ^ A = 0.53GeV, F = (0.30 ± 0.06)GeV^/^ (4.7) 



Combining ( ^^ and ( p.l3|) , the form factors Ai, A2, A3 and V can be calculated as 



functions of T, and sq. In Fig.l, we present our results for the form factors at the zero 
momentum transfer point = 0, which shows the variations of these form factors with 
respect to the Borel parameter T at different threshold energy sq. The T range of interest 
should be similar to that in the light cone sum rule analysis for the B — >■ nlu decay |T^, i.e.. 



T ^ 2.0GeV. It can be seen from Fig.l that, in general consideration, the good stability of 
the form factors exists at the threshold Sq = 2.1 ± 0.6GeV. With such a threshold energy, 
we are now in the stage to evaluate the dependence of the form factors. It should be 
intuitive and convenient to represent these form factors in an algebraic representation. We 
parametrize each form factor in terms of a set of three parameters as follows, 

^(5') = 1 2, f °l r 2/ 2,2 ' (4-8) 

where -F(g^) can be any one of Ai(g^), A2{q^), ^3(5^) and V{q'^). The parameters F{0), ap 
and bp can be fitted from the the sum rule results ( |3.13| ). The results at Sq = 2.1GeV are 
presented in table 1. The form factors as functions of the momentum transfer are also 
shown in Fig. 2. 





F(0) 


ap 


bp 


A, 


0.257 


0.352 


-0.239 


A2 


0.253 


1.090 


0.202 


A3 


-0.253 


1.090 


0.202 


V 


0.134 


1.027 


-0.223 



Table 1. Results of the three parameter fit ( |4.8| ) for the B ^ plv decay 
form factors. These data are fitted from the sum rules ( |3.13| ) at the 
Borel parameter T = 2.0 GeV and the threshold energy sq = 2.1 GeV. 



When the lepton masses are neglected, the differential decay width of i? ^ plv with 
respect to the momentum transfer q^ is 

j-p ri2 IT/ 12 
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with the hehcity amphtudes 



H± = {mB + mp)Ai{q^) T 



Ai/2 



ruB + rrii 



Ho 



2m 



{(m| - ml- q^){mB + mp)Ai{q^) 



X 



rriB + mp 



and 



\ / 2 I 2 2\2 ^22 

A = (m^ + mp — q ) — AmBmp. 



(4.10) 



(4.11) 



So, with the meson masses = 5.28GeV, 
transfer g, 



2 

max 



m 



B + mp — 2mBmp, the integrated width of B 



0.77GeV and the maximum momentum 
plu turns out to be 



T{B plu) 



(10.6±4.0)|Kfe|V 



(4.12) 



The error in eg. ( [4. 121 ) results from the variation of the threshold energy sq = 1-5 — 2.7GeV. 

On the other hand, the branching fraction of p~l^v is measured by CLEO Col- 

laboration pi, Br(5° ^ P'l^T^) = (2.57 ±0.29l[5:i ±0.41) x 10"^ With the world average 



of the lifetime | 



1.56 ± 0.06 ps, one has 
r(S° ^ p-tu) = (1.65 ± 0.80) X 10-^ps-^ 
From eqs. (|4.12|) and (|4.13|) we get 

\Vub\ = (3.9 ±0.6 ±0.5) X 10"^ 



(4.13) 



(4.14) 



where the first (second) error corresponds to the experimental (theoretical) uncertainty. 
Here the theoretical uncertainty is mainly considered from the threshold effects. In gen- 
eral, both the higher twist distribution functions and the QCD radiative corrections may 
modify the parameters and the sum rule results in eg. ( |3.13| ). The modification from higher 
distribution functions is hard to calculate at present because little is known about the high 
distribution functions themselves. It is found from the sum rule analysis that the two-loop 
QCD perturbative corrections may enlarge the constant F by about 25%, and increase A 
at the same time [TH]. To be consistent, one should also include the QCD corrections for 



the correlation function at the same order. As such higher order QCD corrections to the 
correlator have not been considered here, we should take the two-loop QCD corrections of 
the decay constant F as the theoretical uncertainty. Including this uncertainty, our final 
result for \Vub\ is 



\Vub\ = (3.9 ±0.6 ±0.7) X 10" 



This value may be compared with the one obtained in Ref. |12| from B 
using the same approach, 

iKfel = (3.4 ±0.5 ±0.5) X 10~^ 



(4.15) 
ttIu decay by 



(4.16) 



The coincidence of ( [4.15| ) and ( [4.16| ) within their errors proves the consistency of our light 
cone sum rule calculations of heavy to light semileptonic B decays in the framework of the 
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effective theory of heavy quark. We also notice that the theoretical error in ( [4 .151) is larger 
than that in ( 4.16| ), which is not out of expectation since we only include the leading twist 
distribution functions in the B plv calculation. This reflects the importance of a more 
complete study on the p distribution functions. 

The estimate in eq.( |4.15| ) is in agreement with that derived from full QCD calculation 
0, and it is also close to the combined result from the analyses based on different models 
and treatments on B ^ 'k{p)Iv transitions, 

iKfel = (3.25 ± 0.14;°:2i ± 0.55) X 10"^ (4.17) 

which is given by CLEO 



V. SUMMARY 

We have studied B —* plu decay within the framework of effective theory of heavy quark. 
In the effective theory, the relevant matrix elements can be expanded in powers of the inverse 
of the heavy quark mass. At the leading order approximation, the form factors concerned 
in this decay are related to four universal wave functions, which are independent of the 
heavy quark mass mg. Though the HQS loses some predictive power in the heavy to light 
decays, it would be helpful for relating different heavy to light decay channels. For example, 
B ^ plv and D —>■ plu decays are characterized at the leading order of l/m,Q by the same 
set of wave functions, Lj(i = 1, 2, 3, 4). 

The form factors for B plv have been calculated in the effective theory using the light 
cone sum rule approach. The important CKM matrix element \Vuh\ has been extracted by 
comparing the values of integrated width obtained from sum rule calculations and from the 
experimental measurements. The result is 

iKbl = (3.9±0.6±0.7) X 10~l (5.1) 

This result agrees with both the values extracted from the full QCD calculations and that 
from the B — > ttIv decay by using the same approach within the framework of the effective 
theory of heavy quark. This calculation has further shown the reliability of the heavy 
quark expansion and the predictive power of light cone sum rule approach in studying 
heavy to light exclusive decays. We have used B plv decay as an example for concret 
discussion. However, the method is general and the main formulae in this paper can be used 
to other heavy to light vector semileptonic decays after trival modifications such as simple 
replacement of some parameters. 

In this paper, we have considered only the leading twist 2 distribution functions of 
p meson. Both the higher twist and loop corrections and higher order l/mq contributions 
should be included for a more accurate estimation of \ Vub\- It is noted that higher order l/mq 
corrections may have different forms in the usual heavy quark effective theory and the new 
framework of heavy quark effective field theory |T3 HT^] due to the antiquark contributions. 
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FIGURES 





Fig.l(a-d). Variation of form factors with the Borel parameter T for 
different values of the continuum threshold sq. The dashed, solid and 
dotted curves correspond to sq =1.5, 2.1 and 2.7 GeV respectively. Con- 
sidered here is at the momentum transfer = OGeV^. 
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F%"f(a)' 



Frg^2\bj 




Fig.2(a-d). Form factors obtained at the fixed Borel parameter T = 2.0 
GcV. The dashed, sohd and dotted curves correspond to so =1.5, 2.1 
and 2.7 GeV respectively. 
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